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Abstract 

A systematic analysis is presented of the B s — > K form factor / (q 2 ) in the whole range of 
momentum transfer q 2 , which would be useful to analyzing the future data on B s —>■ K decays 
and extracting \V u b\. With a modified QCD light cone sum rule (LCSR) approach, in which the 
contributions cancel out from the twist 3 wavefunctions of K meson, we investigate in detail the 
behavior of / (q 2 ) at small and intermediate q 2 and the nonperturbative quantity fB*9B* b s k {Jb* 
is the decay constant of B* meson and Qb*b s k the B*B S K strong coupling), whose numerical 
result is used to study q 2 dependence of / (q 2 ) at large q 2 in the single pole approximation. 
Based on these findings, a pole model from the best fit is formulated, which applies to the 
calculation on / (q 2 ) in the whole kinematically accessible range. Also, a comparison is made 
with the standard LCSR predictions. 
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A study on heavy-to-light exclusive processes plays a complementary role in the determi- 
nation of fundamental parameters of the standard model(SM) and in the development of QCD 
theory. At present, an important task in the SM is to extract the CKM parameter | V u b \. 
Recently, a new QCD factorization formula |IJ has been proposed for nonleptonic B decays and 
has been applied to discuss phenomenology of B — > tttt, ttK and ttD. This approach, however, 
is not adequate to extracting | V u b | from the relevant data, for the effects of the long distance 
QCD are anyway difficult to control in exclusive nonleptonic processes. Semileptonic B decays 
into a light meson, induced by b — > u transition, are regarded as the most promising processes 
suitable for such a purpose. Nevertheless, in this case precision extraction of | V u b | requires a 
rigorous estimate of the relevant hadronic matrix elements. It is a great challenge, because of 
our inability to deal with nonperturbative QCD effects from the first principle. Heavy quark 
symmetry is less predictive, and also lattice QCD calculation is restricted to a certain kine- 
matical region, for heavy-to-light decays. Although QCD sum rule method has been being an 
effective QCD-based approach to nonperturbative dynamics, with its extensive uses in phe- 
nomenology, the resulting form factors behave very badly in the heavy quark limit rriQ — > oo, 
the reason being that one omits the effects of the finite correlation length between the quarks in 
the physical vacuum. In order to overcome the disadvantage, an alternative to the "classical" 
sum rule method, QCD light-cone sum rules (LCSR) , has been developed and has proven to 
be an advanced tool to deal with heavy-to-light transitions. There exist a lot of its applications 
in the literature. For a detailed description of this method, see J|. However, a problem to need 
solving is how to control the pollution by higher twist, especially twist 3 wavef unctions, which 
are known very poorly and whose influence on the sum rules is considerable in most cases. In 
||, an improved LCSR approach has been worked out, to eliminate twist 3 wavefunctions 
and enhance the reliability of sum rule calculations, and has been applied to reexamine heavy- 
to-light form factors in the region of momentum transfer < q 2 < m| — 2mi ) A.QCDi where the 
operator product expansion (OPE) goes effectively in powers of small light-cone distance x 2 . 

Most of previous works 0, ||, @, 0, @| are devoted to discussing B — > ir, p semileptonic 
transitions within the context of LCSR, with the aim to extract | V u b |. A study on B s — > K 
semileptonic processes is equally important. As compared with the case B — > tt, p, however, 
the B s — > K form factors are more difficult to evaluate, for SU(3) breaking corrections to the 
twist 3 wave functions of K meson have not been investigated completely in the literature. 
Explicitly, this problem can be avoided in our approach |4|, [5| . On the other hand, to calculate 
the semileptonic widths one must find another way to estimate the form factors at the large 
momentum transfer m 2 — Ivti^Kqcd < <? 2 < ( m B s — m K) 2 ■ in the letter, we investigate the 
B s —>■ K form factor f(q 2 ) at the total momentum transfer with the improved LCSR and a 
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pole model. 

Let us start with the following definition of the B s — > K form factors f{q 2 ) and f(q 2 )'- 

(K{p)\uj,b\B(p + q)) = 2f(q 2 )p, + f(q 2 )q,. (1) 

For B s — > KiUf transitions, as I = e, p we can neglect the contributions from f(q 2 ) due to 
the smallness of m e ^ and therefore only the form factor f(q 2 ) is relevant. It can precisely be 
represented as 



f{q2) = Jb* 9b* BsK + f P(a)da 



2tob.(1 -g 2 /mL) J l — q 2 /a 

2\ , 17 /„2\ 



= F G (q 2 ) + F H (q 2 ), (2) 

with f B * being the decay constant of B* meson, mg. the B* meson mass, Qb*b s k the strong 
coupling defined by 

(B*(q,e)K(p)\B s (p + q)) = -g B * BsK {p ■ e), (3) 

and p(a) a spectral function with the threshold a^. Obviously, Fc{q 2 ) stands for the contribu- 
tion from the ground state B* meson, which describes the principal behavior of f(q 2 ) around 
q 2 = q^ax, and F B (q 2 ) parametrizes the higher state effects in the B* channel. As we have 
known, the form factor f(q 2 ) may be estimated for the small and intermediate momentum 
transfers by means of LCSR, and also the nonperturbative parameter f B *g B * Ba K is accessible 
within the same framework. Accordingly, modelling the higher state contributions by a certain 
assumption and then fitting (2) to its LCSR result /lc(? 2 ) m the region accessible to the light 
cone OPE, we might derive the form factor f(q 2 ) in the total kinematical range to a better 
accuracy. For this purpose, we follow the procedure in [|], [5| and consider a chiral current 
correlator used for a LCSR sum rule calculation on /lc(<? 2 ) and j b*Qb*b s Ki 

n M (p,g) = i Jd 4 xe^(K(p)\T{u(xh,(l + 75 )6(x), 6(0)i(l + 75 ) s (0)}|0> 

= F(q 2 ,(p + q) 2 )p, + F(q 2 ,(p + q) 2 ) qfM (4) 

Inserting complete sets of the relevant intermediate states \B H ) in (4) and using the definition 
(0|s«7 5 6|5 s ) = m ^ f Bs and {0\u^^b\B*) = m B *f B *e^, we have the two hadronic representa- 
tions of the invariant function F(q 2 , (p + q) 2 ), 

,h, 2 / , \ 2 \ 2fLc(q 2 )m 2 Bs fB s 7 Pi (s) 



F 1 (q,( P + q) )" fa + ma )( m % -(p'+q)*) +J s _ {p + q) 2 ds ( 5 ) 



3 



,// 2 m B s m B* Ib 3 Ib* 9b*B s K 



W,(p + ?r) 



f f P 2 [si, s 2 )V{ Sl - s )V{s 2 - s ) 
+ J J (s,-a*)(s*-(v + a)*) dSldS2 - (6) 



(m b + m s )(m%, - q 2 ){m 2 Bs - (p + q) 2 ) 
p 2 (si,s 2 )Q(s 1 - s' )Q(s 2 - s ) 
(si - q 2 ){s 2 - (p + q) 2 ) 

Several definite interpretations for (5) and (6) are in order. The two dispersion integrals include, 
in addition to the contributions of the resonances carrying the same quantum numbers as the 
corresponding ground states in the pole terms, the effects due to the relevant orbit-excited 
B mesons. Taking it into account that these orbit-excited states are far from the lowest B s 
and B* mesons, and the lowest two of them are slightly below the first excited B s and B* 
mesons in mass, their contributions can effectively absorbed into a dispersion integral so that 
thresholds s and s' should correspond to the squared masses of the lowest + B s and 1 + 
B mesons respectively. On the other hand, the vector current u^^b and axial-vector current 
u'y^'y^b couple also to + and 0~ B mesons, respectively, which should be considered in (6). The 
invariant function, however, does not receive such a contribution as we have checked. Therefore, 
it is safe to separate the hadronic expression F^(q 2 , (p + q) 2 ) into a pole term and a dispersion 
integral. 

The task left is to calculate the correlator in QCD theory in order to obtain the desired sum 
rules. To this end, we work in the large space-like momentum regions: (p + q) 2 <^ for the 
ficiq 2 ) case and q 2 <C 0, (p + q) 2 for the fB s fB*9B* b s k case, so that the light cone OPE 
can be used for the correlator under consideration. After contracting the b quark operators, we 
encounter some nonlocal matrix elements, which can systematically be expanded in powers of 
the deviation from the light cone x 2 = by defining the relevant light cone wavefunctions of K 
meson classified in terms of twist. The explicit parametrizations of all those can be found in 
the literature, and here are no more given for simplicity. A long but straightforward calculation 
yields the light-cone QCD form of the invariant function F(q 2 , (p + q) 2 ), 



F QCD {q 2 ,{p + q) 2 ) = 2f K m b < 



l 

du 



o 



<p K {u) Bml\gi(u) - G 2 (u)\ 2ug 2 {u) 



m b ~ (q + u p) 2 [ m b ~ (q + up) 2 } 3 



m l ~ (q + u pY 



i 



' [m( - {q + (3p) 2 \ 2 1 



o 



to twist 4 accuracy. Here <Pk(u) is the twist 2 wavefunction, while the others have twist 4; 
the parameter (5 — a± + aa^ and Don = da\da 2 daz8 (1 — cki — a 2 — a 3 ) . At this point, we 
put once again an emphasis on that differently from the existing LCSR calculations, the twist 
3 wavefunctions make precisely a vanishing contribution to the correlator we choose. This is 
essentially important to enhancing precision of the LCSR calculation. 



As usual, we need to make the Borel improvements on the theoretical expression: F® CD (q 2 , (p+ 
q f) _^ Fi QCD (q 2 , M 2 ), F QCD {q 2 , (p + q) 2 ) -> F 2 QCD (M 2 , Mf ), and then match them onto the 
individual Borel improved hadronic forms. Invoking the quark-hadron duality ansatz, the final 
sum rules for fic(q 2 ) and gs* b s k read respectively: 



f(<l 2 



Tnb( m b + m s 



-fxe 



m B s 
M' 2 



1 du 



A U 



m^-(q 2 -um 2 K )(l-u) 



[<Pk ( u ) 



Ami 
u 2 M 4 



uM 2 



g 2 (v)dv 1 + 



o 



m 2 b + q 2 " 
uM 2 



e(3-A) ™i-(i 2 -e™ 2 K )^-0) 

+ da Dai 1 : ^ e ^ 



„ 2 — ^£ 

— 4mie a?2 



— 2e a? 3 



/? 2 M 2 
1 



q 



2\2 



1 + 



+ q 2 



(s - q 2 ){m 2 b - q T 



sp-q 

M 2 
92(A) - 



<?i(A) 



1 



1 



dgi(A) 



(sq — q 2 ){ml — q 2 ) du 



(m 2 b - q 2 ) 



1 + 



ml + q 2 
ml — q 2 



1 + 



go - q 

M 2 



g 2 (v)dv 



Jb s Jb* gB*B s K 



2 1r 

2m b (m b + m s )f K m B s ^B* ( 2 

2M ' 1 • 



M 



m%m B * 



2 1 2 



— e a? 



¥>*(l/2) 



2 12 

-e m 2 



, (1 / 2 )-f 



1/2 



^(1/2) - / </ 2 (u)du 



1/2 

+ / <iai 



l-a 



da? 



where A 



and M 



l/2-Qi 
M 2 M 2 



Q'3 



2^_ L (o! i ) + 2p ± (ai) - (p\\(ai) - £||(ar f ) 



. In the derivation of (9), we have taken M 2 



(9) 



M? due 



to the fact that B s and B* mesons are nearly degenerate in mass, which renders the continuum 



subtraction reduce to a simple replacement e 
2 term. 



go 

■ e m 2 for the leading twist 



Turning now to the numerical discussions on the sum rules. B channel parameters entering 
the sum rules are the b quark mass m&, B meson masses wb s and mg., decay constants fs s 
and fs*, and threshold parameter so- We take ms, = 5.369 GeV, mg. = 5.325 GeV and 
wife = 4.8 GeV. As for the decay constants /# a and fs*, we have to reanalyze them in the two- 
point QCD sum rule approaches 0, |5| with chiral current correlators, to keep a consistency with 
the sum rules in question. The results are found to be fg s = 0.142 GeV and = 0.132 GeV, 
as the threshold parameter sq = 34 GeV 2 corresponding to the mean value of squared masses 
of the lowest + B s and 1 + B mesons. For the decay constant of K meson and mass of s 
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quark, we use fx = 0.16 GeV and m s = 0.15 GeV. The important point is to specify the 
set of the light-cone wave functions of K meson. Unlike the case of ir meson, SU(3) breaking 
effects need considering for the distribution amplitudes of K meson. In the work, we use the 
model presented in |J for the leading twist wavefunction, which is based on an expansion over 
orthogonal Gegenbauer polynomials with coefficients determined by means of QCD sum rules. 
The explicit expression is 



(Pk(u) = 6u(l — u) \ 1 + 1 



0.5(2w 



1 + 1.2[(2 M -1) 2 -^] 



(10) 



at the scale = yrn 2 B s — m 2 , measuring the mean virtuality of the h quark. For the twist-4 
wave functions, we neglect the SU(3) breaking effects and utilize the same forms as those of 7r 
meson investigated in JTT)| |. 

Having fixed the input parameters, one must look for a reliable range of the Borel parameters 
M 2 and M , which can be determined by the standard procedure. The fiducial intervals are 
found to be 8 GeV 2 < M 2 < 17 GeV 2 , depending slightly on q 2 , for q 2 = - 17 GeV 2 and 
5 GeV 2 < M < 10 GeV 2 . In the two "windows", the twist 4 wavefunctions contribute less 
than 9% and 7%, and the continuum states at the levels lower than 25% and 22%, respectively. 
The sum rule results for /lc(<? 2 ) show a weak dependence on M 2 up to q 2 = 17 GeV 2 , varying 
between ±3%— ±5% relative to their central values. For the product fB a fB*9B* b s k, the resulting 
sum rule is f B s f b* 9 b* b s k = 0.55 GeV 2 , the uncertainty due to M being ±4%. Taking its 
central value, we get gs*B a K = 29. To evaluate better the B pole contribution in (2), however, 
we would give a direct sum rule result for /b*9b* b s Ki which can be obtained utilizing the 
analytic form instead of the numerical result for the two-point sum rule for f^ s in (9). The 
result is fB*9B* b s k = 3.57 — 4.19 GeV, depending on the Borel parameters. The sum rule 
prediction /lc(9 2 )> together with that from the B* pole approximation, is illustrated in Fig. 1. 
It is explicitly demonstrated that a perfect match between them appears at q 2 ~ 15 — 20 GeV 2 . 

The influence on the sum rules should be investigated in detail from several important 
sources of uncertainty: the twist 2 distribution amplitude ip K {u), b quark mass m h , decay 
constants fg a and fg*, and threshold parameter So- Concerning the light cone wavefunction 
(Pk{u), there are some determinations other than that in (10) in the literature. To investigate 
the sensitivity of the sum rules to the choice of the non-asymptotic coefficients in <^(w), we 
consider the two models suggested in and ||, and confront the resulting sum rules with our 

ones. If adopting ipxiu) in 0, the resulting changes amount to —8% 9% for the fLc{l 2 ) 

case and to ±5% for the Jb*9b* b s k case. The almost same situation exists for that used in 0. 
Therefore, the uncertainties caused by <Pk{ u ) ma y be estimated at a considerably small level. 



As for the B channel parameters m^, fs a , Jb* and Sq, considering a correlated variation in the 
individually allowed ranges would give sufficient information on the uncertainties induced by 
them. This can be done in such a way where letting mj vary from 4.7 to 4.9 GeV, we observe 
the behavior of /lc(<? 2 ) and /b*9b* b 3 k by requiring that the relevant decay constants take only 
the best fitting values. We find that such an effect amounts to 6% and 5%, respectively. At 
present, the total uncertainties in ficil 2 ) an d f b* 9b* b s k can respectively be estimated to be 
20% and 18%, by adding linearly up all the considered errors. 

It is important and interesting to make a comparison of our sum rule results and those 
from the standard LCSR based on the corrector of vector and pseudoscalar currents, which are 
easy to obtain using the twist-3 wavefunctions suggested in , leaving the twist-4 distribution 
amplitudes unchanged and making a corresponding replacement of the other relevant input 
parameters in (79) and (44) of the second reference in ||. We observe that the standard 
approach gives the same matching range as in our case and the resulting deviations from our 

predictions turn out to be between —10% 15%, depending on q 2 , in the total kinematically 

accessible region. This denotes that both approaches are essentially compatible with each other 
within the available errors. 

With the yielded findings we would give a specific parametrization for f(q 2 ) applicable to 
the whole kinematical region, which is helpful for the future practical application. Assuming 
the higher state contribution in (2) to obey Fn{q 2 ) = a/ (1 — bq 2 /m 2 B * — cq A /m Bt ), we have a 
pole model for f(q 2 ), 

9b*b s k a ^ 

2mg» (1 — q 2 /m 2 B ,) 1 — bq 2 /m 2 B * — cq 4 /m Bt 

The parameter a can easily be fixed at —0.07, using the central values of /lc(0) and }'b*9b*b s k- 
In the region q 2 = - 18 GeV 2 , the best fit of (11) to hc{q 2 ) yields b = 1.11, and c = -8.33. 
The resulting q 2 dependence of f(q 2 ) is demonstrated in Fig.l too, for a comparison. It turns 
out that the fitting results reproduce precisely the LCSR prediction up to q 2 = 18 GeV 2 and 
support considerably the single pole description of the B s — > K form factor f(q 2 ) at large q 2 . 

Also, it is worthwhile to look roughly into SU(3) breaking effects in heavy-to-light decays 
by considering the ratio of the derived B s — > K form factor over the corresponding B — ► ix one. 
The B — ► 7i form factor has already been obtained for small and intermediate q 2 in the improved 
LCSR approach in Q. Using all the same method as in present case, we can understand its 
behavior at large q 2 and further get a parametrization holding for the total kinematical range. 
For the common kinematical region to the two processes, the resulting ratios, a comparable 
result 1.05 — 1.15 with that from the standard approach, favor a small SU(3) breaking effect. 

We have given a detailed discussion on the B s — > K form factor f(q 2 ) in the whole kine- 



matical region. To avoid the contamination with the twist 3 wavef unctions, in which SU(3) 
breaking corrections have not been analyzed systematically, an improved LCSR approach with 
some kind of chiral current correlator has been applied to estimate the form factor /lc(? 2 ) a t 
small and intermediate q 2 . The nonperturbative quantity j b*Qb*b 3 k-, an important input in the 
B* pole model for f(q 2 ), has also been calculated within the same framework and the sum rule 
result has been adopted to study the behavior of f(q 2 ) at large q 2 . We find that the resulting 
fhcil 2 ) matches quite well with the estimate from the B* pole model at q 2 = 15 — 20 GeV 2 . 
A comparison shows that our predictions are in basic agreement with those from the standard 
LCSR. Based on our findings, a pole model for f(q 2 ) has been worked out, which is applicable 
to the total kinematically accessible region. The results presented here would be used as ana- 
lyzing the future data on B s — > K decays and extracting | V uh |. A future lattice calculation 
of the B s — > K form factors, which is available for large q 2 , will provide a direct test of our 
predictions. The same approach applies also to discuss other heavy-to-light processes. 
This work is in part supported by the National Science Foundation of China. 
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Figure 1: The B s K form factor f(q 2 ) in the total kinematical range. The solid line denotes 
the LCSR result /lc(<? 2 )> which is reliable for < q 2 < 17 GeV 2 . The dotted line expresses 
the B* pole prediction suitable for large q 2 . The best fit of Eq. (11) to fhc{q 2 ) is illustrated 
by the dashed line. It should be understood that the plotted curves correspond to the central 
values of all the relevant parameters. 
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